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We extend a model of four-dimensional simplicial quantum gravity to in- 
clude degenerate triangulations in addition to combinatorial triangulations tra- 
ditionally used. Relaxing the constraint that every 4-simplex is uniquely defined 
by a set of five distinct vertexes, we allow triangulations containing multiply 
connected simplexes and distinct simplexes defined by the same set of vertexes. 
We demonstrate numerically that including degenerated triangulations substan- 
tially reduces the finite-size effects in the model. In particular, we provide a 
strong numerical evidence for an exponential bound on the entropic growth of 
the ensemble of degenerate triangulations, and show that a discontinuous crum- 
pling transition is already observed on triangulations of volume N4 rj 4000. 

PACS numbers: 04.60.-m,Nc 05.70.Fh 02.70.Lq 

Discretized models of four-dimensional Euclidean quantum gravity, known as sim- 
plicial gravity or dynamical triangulations, have received ample attention in recent 
years, the hope being that in a suitable scaling limit they might provide a sensible 
non-perturbative definition of quantum gravity. In the simplicial gravity approach 
the integration over metrics is replaced by a sum over an ensemble of triangulations 
constructed by all possible gluings of equilateral 4-simplexes into closed (piece-wise 
linear) simplicial manifolds (see e.g. Ref. 0, §])■ The regularized Euclidean Einstein- 
Hilbert action is particularly simple; it can be taken to depend on only two coupling 
constants, fj, and k, related to the cosmological and the inverse Newton's constants. 
The coupling constants are conjugate to the volume — the number of 4-simplexes 
— and the number of triangles in a given triangulation respectively. The regularized 
grand-canonical partition function thus becomes: 



The sum is over all distinct triangulations T € T , iVj is the number of i-simplexes in 
a triangulation T and Ct denotes its symmetry factor — the number of equivalent 
labeling of the vertexes. 
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Extensive numerical simulations have established that the model Eq. (||) has a 
strong-coupling (small k) crumpled phase and a weak-coupling (large k) elongated 
phase, separated by a discontinuous phase transition. In the crumpled phase the 
internal geometry collapses and is dominated by a novel singular structure — two 
singular vertexes connected to an extensive fraction of the total volume and joined 
by a sub-singular edge Q. The elongated phase, on the other hand, is dominated 
by branched polymer like triangulations, i.e. bubbles glued together via small necks 
into a tree-like structure. 

In Eq. (jl]), T denotes a suitable ensemble of triangulations included in the par- 
tition function. Different ensembles are defined by imposing various restrictions on 
how the simplexes are glued together. Provided this leads to a well-defined parti- 
tion function, and as long as the difference is only at the level of discretization, one 
expects different choices of T to lead to the same continuum theory in the thermody- 
namic limit. This is known to be true in two dimensions where models of simplicial 
gravity corresponding to different choices of T are soluble as matrix models ||. 
Even in two dimensions, however, for the partition function to be convergent the 
topology of the triangulations included in the sum Eq. (|l|) must be fixed, regardless 
of the ensemble used. As the same most likely is true in higher dimensions, in this 
letter we only consider triangulations of fixed spherical topology. 

All simulations of four-dimensional simplicial gravity have, as of yet, used an 
ensemble of combinatorial triangulations Tq. In a combinatorial triangulation every 
D-simplex is uniquely defined by a set of (D + 1) distinct vertexes — it is said to 
be combinatorially unique. In this letter we study a larger ensemble of degenerate 
triangulations where we relax this constraint and allow distinct simplexes to be 
defined by the same set of vertexes. This includes two simplexes with more than 
one face in common. We do, however, retain the restriction that every 4-simplex is 
defined by a set of five distinct vertexes, i.e. we exclude degenerate simplexes. Clearly 
Tc C 7b- This corresponds to an ensemble of restricted degenerate triangulations 
as defined in Ref. [||. 

The benefits of using a larger ensemble of triangulations are well known from 
simulations of two-dimensional simplicial gravity which have demonstrated that less 
restricted the triangulations are translates into smaller finite-size effects |?|]. Re- 
cently the same observation has been made in three dimensions |J. As simulations 
of four-dimensional simplicial gravity are notoriously time-consuming, primarily due 
to the large volumes needed to observe any "true" infinite-volume behavior, any re- 
duction in the finite-size effects is of great practical importance. 

In this letter we show that including degenerates triangulations in simulations 
of four-dimensional simplicial gravity likewise leads to reduced finite-size effects. 
This reduction is most pronounced in the crumpled phase of the model where the 
free-energy density of the canonical (fixed volume) ensemble — the pseudo-critical 
cosmological constant fi c (N^) — converges very rapidly to an infinite- volume value. 
This in turn implies an exponential bound on the entropic growth of the ensemble of 
degenerate triangulations, something that has been the subject of some controversy 
in the past for combinatorial triangulations [||, [| . Although we observe qualitatively 
the same phase structure as with the model Eq. (||) restricted to combinatorial tri- 
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angulations, there are some dissimilarities. A discontinuous phase transition, sepa- 
rating the elongated and the crumpled phases, is already observed on triangulations 
of relatively modest size, N4 ~ 4000, compared to combinatorial triangulations were 
a volume of N4 ~ 32.000 is needed. And while the crumpled phase is still dominated 
by a singular structure, for degenerate triangulations this corresponds to a gas of 
sub-singular vertexes rather than to only two singular vertexes. 

We have simulated the model Eq. (|l]) using degenerate triangulations on volumes 
up to 32.000 4-simplexes using Monte Carlo methods. As customary we work in a 
quasi-canonical ensemble of spherical manifolds with almost fixed N4: 

Z(p,n;N 4 ) = ^e-^-^-^Q^), (2) 

where VLn 4 (k) = Stst^) ex P( K -^2) ^ s ^ e canonical partition function. As there 
do not exist ergodic volume conserving local moves, hence the canonical ensemble 
cannot be simulated directly, we must allow the volume to fluctuate. The quadratic 
potential term added to the action ensures, for an appropriate choice of 5, that these 
fluctuations are small. 

In the simulations the triangulation space is explored using a set of local geometric 
changes, the (p, g)-moves. In a (p, q)-move, where p = D + l — q, a (g— l)-subsimplex 
in the triangulation is replaced by its "dual" (p— l)~subsimplex. For combinatorial 
triangulations the (p, q)-moves are known to be ergodic for D < 4 [j[C|| . To demon- 
strate that the same holds true for degenerate triangulations we observe that, just 
as in three dimensions ||, every set of combinatorially equivalent simplexes, or 
sub-simplexes, can be made distinct by a finite sequence of the (p, g)-moves. Thus 
every degenerate triangulations can be reduced to a combinatorial one. In addi- 
tion the local nature of the (p, (/)-moves prohibits the creation of pseudo-manifolds 
in the simulations, i.e. triangulations containing vertexes with a neighborhood not 
homeomorphic to the D-ball. 

From a practical point of view simulating degenerate triangulations is actually 
simpler than simulating their combinatorial counterpart as one avoids the non-local 
manifold checks necessary to exclude combinatorially equivalent simplexes. For 
combinatorial triangulations these checks are the most time-consuming part of the 
simulations |Oj. This simplification of is particularly beneficial in the crumpled 
phase where the singular structure dominates; in this phase we observe a tenfold 
reduction in the effective auto-correlation times (measured in "real" time) when 
using degenerate instead of combinatorial triangulations. In the branched polymer 
phase, on the other hand, the auto-correlation times appear comparable for the two 
ensembles. 

The real benefit of using degenerate triangulations is the reduction of geometric 
finite-size effects. This reduction is most striking for the volume dependence of the 
pseudo-critical cosmological constant, p^(N±), which we shown in Figure 1. For 
comparison we also include the corresponding values, p^? (N4), for combinatorial 
triangulations. For degenerate triangulations we observe a rapid convergence to an 
infinite volume value p. This can be quantified by comparing the fit of pf? (N4) to 
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Figure 1: The pseudo-critical cosmological constant n®(N±) for an ensemble of degenerate 
triangulations, together with fits assuming a power-law convergence, Eq. (|^) (solid curve), 
or a logarithmic divergence, Eq. (|J) (dashed curve). Also included are the corresponding 
values, /ig (.ZV4), for an ensemble of combinatorial triangulations. 

Table 1: The parameters in the fit of the pseudo-critical cosmological constant /Lt^ (-/V4) to 
Eqs. (||) and (||) respectively. Measurements on volumes N 4 — 400 to 25.600 are included 
in the fits. 



\i 7 X A /d.o.f. 


Eq. (1 
Eq. ([ 


j) 2.556(3) 0.55(5) [3.8] 
j) 2.385(4) [117] 



two different functional forms: a weak power-law convergence, 

McW = A + ^7> (3) 
4 

or a logarithmic divergence, 

Mc(AT 4 ) = /2 + 6 , logJV 4 . (4) 

The fit parameters and the quality of the fits are shown in Table 1. In contrast to 
combinatorial triangulations for degenerate triangulations there is no comparison in 
the quality of the fits; the latter, which corresponds to a divergent partition function 
Eq. (|l|), is ruled out by a x 2 /(d.o.f.) f» 117. For combinatorial triangulations, on 
the other hand, it is difficult to use the quality of the fit of n^{N^) to either Eq. @ 
or Eq. ([|) to distinguish between those two scenarios (see e.g. Ref. ||). 
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Figure 2: The (normalized) distributions of local volumes, P(q), for degenerate triangula- 
tions. This is for 400 < JV4 < 3200, and for k = 0. (Insert) The number, n s , of sub-singular 
vertexes versus volume. 

The importance of this result is that it provides strong numerical evidence for an 
exponential bound on the entropic growth of the ensemble of degenerate triangula- 
tions as a function of volume — a necessary condition for a well-defined partition 
function Eq ([!]). And, as Tq £ 7b > this implies an exponential bound on the number 
of combinatorial triangulations as well. 

The origin of the large finite-size effects present in simulations with combinatorial 
triangulations lies in the nature of the quantum geometry in the crumpled phase. 
As stated, the partition function is dominated by triangulations characterized by 
two singular vertexes connected by a sub-singular edge. A singular vertex has a 
local volume q — the number of 4-simplexes containing that vertex — which grows 
linearly with the volume of the manifold, while the sub-singular edge has a local 
volume which grows like N£, a ~ 2/3 |J. However, for combinatorial triangulations 
this singular structure only dominates on large enough volumes, on small volumes 
triangulations with only one singular vertex have larger entropy. This results in a 
cross-over behavior in the fractal structure at A4 ~ 1000, as can be observed in 
Figure 1. 

For degenerate triangulation the crumpled phase is likewise dominated by a sin- 
gular structure. This is evident from the probability distribution P(q) of the local 
volumes which contains an isolated peak in the tail. This is shown in Figure 2 for 
k = 0. However, the distribution P(q) differs in two respects from the corresponding 
distribution measured on combinatorial triangulations |J: 

(a) The peak corresponds to not just two singular vertexes but rather to several 
sub-singular vertexes, i.e. vertexes with local volumes that scale like N", a < 1. 
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A rough estimate yields a ~ 0.9. The number of these sub-singular vertexes, n s 
increases logarithmically with the volume as is shown in the insert in Figure 2. 
This suggests that the crumpled phase is dominated by a gas of sub-singular 
vertexes. 

(b) For each volume, P(q) effectively separates into two distinct distributions 
depending on whether the local volume is even or odd. 

It is not clear though how much significance should be attached to this difference 
in the singular structure. Due to the collapsed nature of the internal geometry it is 
unlikely that any sensible continuum limit exists in the crumpled phase, hence there 
is no reason to expect identical scaling behavior for the two different ensembles. 
The details of the discretization may still be important in the thermodynamic limit 
for k < k c . It is, however, worth noticing that for degenerate triangulations we do 
not observe any change in the fractal structure as the volume is increased as for 
combinatorial triangulations. 

Additional evidence of a collapsed intrinsic geometry in the crumpled phase comes 
from the (absence of) volume scaling of the simplex-simplex distribution s(r), i.e. 
the number of simplexes at a geodesic distance r from a marked simplex. Using the 



scaling ansatz s(r) = n\ 1 ^ dli F(x), where x = r/N^ dli [12, [L3|], we tried to collapse 
distributions s(r) measured on different volumes onto a single scaling curve. This 
was though not possible with an "acceptable" collapse (with x 2 /(d.o.f.) of order 
unity); moreover, the estimate of the fractal dimension du appeared to increase 
with the volume. From this we conclude that dn = oo in the crumpled phase. 

We have also investigated the phase structure of the model for non-zero values 
of the inverse Newton's constant k. As for combinatorial triangulations we observe 
a phase transition to an elongated phase at k c ~ 1.7. To establish the nature of 
the phase transition we have studied the Monte Carlo time-series of the energy 
density, no = N0/N4, in the critical region. We show an example of one such time- 
series in Figure 3 together with the corresponding histogram. This is for volume 
A4 = 4000 and n = 1.673. The histogram shows a clear double-peak structure 
characteristic of a discontinuous phase transition. Note that in order to observe the 
corresponding two-state signal using the combinatorial ensemble, triangulations of 
volume N 4 ~ 32.000 are needed @]. 

For k > k c the model is in an elongated or branched polymer phase. This we 
have established by measuring the fractal dimensions du and the spectral dimension 
d s for k = 2. The former is determined from the scaling of the simplex-simplex 
distribution, the latter from the return probability of a walker on the dual graph, 
p(t) ~ t~ d s/ 2 Including measurements on volumes A4 = 400 to 1600, we get 



dn = 1-9(2) and d s = 1.32(5), in excellent agreement with dn = 2 and d s = 4/3 as 
expected for branched polymer. 

In this letter we have demonstrated that including degenerate triangulations in 
simulations of four-dimensional simplicial gravity has many potential advantages 
over the model restricted to combinatorial triangulations. This agrees with the same 
observations previously made in both two and three dimensions. The chief benefit is 
the reduction in geometric finite-size effects manly due to an enlarged ensemble - 
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with a larger triangulation-space the infinite- volume fractal structure is more easily 
approximated on small volumes. The most important result presented in this letter 
is a strong numerical evidence for an exponential bound on the entropic growth of 
the canonical partition function Oat 4 . A more practical result is the observation of 
a discontinuous phase transition on triangulations of relatively modest size. 

A natural extension of the work presented in this letter is to investigate how the 
phase structure is affected if the model Eq. ([l]) is changed either by using a modified 



measure |14| or by adding matter fields [15]. For combinatorial triangulations it has 
recently been observed that this can substantially alter the phase structure and, for 
a suitable modification, a new crinkled phase appears Ipf] . If this observed phase 
structure corresponds to a genuine change in the continuum behavior of the model 
Eq (||), one expects on basis of universality that it should be independent of the 
ensemble of triangulations used. 
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